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We study a chaotic quantum transport in the presence of a weak spin-orbit interaction. Our 
theory covers the whole symmetry crossover regime between time-reversal invariant systems with 
and without a spin-orbit interaction. This situation is experimentally realizable when the spin-orbit 
interaction is controlled in a conductor by applying an electric field. We utilize a semiclassical 
approach which has recently been developed. In this approach, the non-Abelian nature of the 
spin diffusion along a classical trajectory plays a crucial role. New analytical expressions with one 
crossover parameter are semiclassically derived for the average conductance, conductance variance 
and shot noise. Moreover numerical results on a random matrix model describing the crossover 
from the GOE (Gaussian Orthogonal Ensemble) to the GSE (Gaussian Symplectic Ensemble) are 
compared with the semiclassical expressions. 

PACS numbers: 73.23.-b, 05.45.Mt, 03.65.Sq, 05.40.-a 



I. INTRODUCTION 

A chaotic quantum transport of an electron in a cav- 
ity is caused by either implanted impurities or bumpy 
boundaries, and provides directly measurable quantum 
signatures of chaoa^, such as the conductance variance. 
Universal aspects of a chaotic transport have been inves- 
tigated by means of the random matrix theory (RMT) 2 . 
In the RMT, quantum systems are classified into sym- 
metry classes. A chaotic system with time-reversal sym- 
metry is described by the Gaussian orthogonal ensemble 
(GOE). When the time- reversal symmetry is broken by 
applying a magnetic field, the Gaussian unitary ensemble 
(GUE) becomes a suitable model. If a system with time- 
reversal symmetry has a spin-orbit interaction, one needs 
to employ the Gaussian symplectic ensemble (GSE). 

We consider the case that two leads are attached to a 
cavity and the number of the lead channels are Ni and 
N 2 . An electron transport in the cavity is described by 
the scattering matrix 3 -^. Replacing the scattering matrix 
by a random matrix, the RMT phenomenologically pre- 
dicts the average conductance G, conductance variance 
VarG, and shot noise P at zero temperature as££ 
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with ^ = ^1+^2. Here j3 — 1,2, and 4 correspond to 
the GOE, GUE, and GSE symmetry classes, respectively. 
These expressions include the contributions from the spin 
degrees of freedom 7 . The constants Go and Po are Go = 
e 2 /(Trh) and Po = 2e 3 |U|/(7r?i), respectively, where e is 
the unit electric charge and V is the bias voltage. If N± is 
equal to N 2 and very large, the leading term of the shot 
noise is insensitive to a change of the symmetry. 

When a very weak magnetic field is applied to the cav- 
ity, the time-reversal symmetry is only partially broken. 
In this case, a crossover from the GOE to GUE is ob- 
served. This GOE-GUE crossover regime can also be an- 
alyzed by a parametric RMT model, and analytic predic- 
tions describing a chaotic quantum transport are known 8 . 
Recently, a chaotic transport in the GSE-GUE crossover 
regime was also studied within the RMT framework^. In 
this regime, a very weak magnetic field breaks the time- 
reversal symmetry of a system with a spin-orbit interac- 
tion. 

The aim of this paper is to study another case, the 
crossover from the GOE to GSE, in which the system has 
a very weak spin-orbit interaction preserving the time 
reversal symmetry. In the experimental point of view, 
the GOE-GSE crossover can be realized, if the spin-orbit 
interaction (or Aharonov-Casher effect) is controlled by 
applying an electric field in a chaotic conducto r 10 ' 11 . In 
this case, a parametric RMT model is also know n 12 ' 13 , 
and the diagrammatic perturbation theory has been 
used to evaluate some transport properties^—. Here 
we employ a semiclassical approach which has recently 
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been developed^—. In a semiclassical evaluation, the 
transmission amplitude is treated by the path-integral 
method, where all the classical paths must in principle 
be taken into account. However recent studies clarified 
that almost the same but partially time reversed pairs of 
classical trajectories contributed to the conductance^— , 
so that the calculation was greatly simplified. Then it 
was shown that the semiclassical approach could pre- 
cisely reproduce the RMT predictions (JTJ-lJSJSl-^. More- 
over, when a similar approach is applied to the paramet- 
ric spectral correlations in the GOE-GUE, GUE-GUE, 
GOE-GOE and GSE-GSE regimes, it can also reproduce 
the RMT predictions^ 3 ^. 

Thus the semiclassical approach has become a practi- 
cal tool to find a new prediction, even if the RMT anal- 
ysis is difficult. As this approach was already applied 
to the parametric spectral correlations in the GOE-GSE 
crossover regime 3 ^*, we naturally expect that it can be 
used in the analysis of a transport. 

Considering the non-Abelian nature of the spin diffu- 
sion along the classical trajectories, we extend the semi- 
classical technique to derive analytic expressions for the 
transport properties. Our results on the average conduc- 
tance, conductance variance and shot noise are given in 
eqs. (HU), (EU), and ([35]). The crossover from the GOE 
to GSE is controlled by one parameter depending on the 
diffusion constant of the spin. The GOE and GSE results 
are reproduced in the limiting cases of the parameter. 

This paper is organized in the following way. In Sec. II, 
a semiclassical expression of the transmission amplitude 
is presented. We put a stress on the statistical aspects 
of the expression. In Sec. Ill, using the semiclassical ex- 
pression, we calculate the average conductance, conduc- 
tance variance and shot noise. In Sec. IV, these results 
are compared with numerical calculations on a random 
matrix model. We finally summarize the paper in Sec. 
V. 



II. SEMICLASSICAL EXPRESSION OF THE 
TRANSMISSION AMPLITUDE 

The semiclassical theory employs the transmission am- 
plitude t ai ^ a21 which represents the propagator of a wave 
packet from the channel a\ in one lead to a 2 in another 
lead. Bolte and Keppeler derived a semiclassical expres- 
sion of the transmission amplitude with spin variables 32 : 

ta u a 2 ~ J2 A a ^ a e^'\ (4) 

a:ai — >a-2 

where Th is the Heisenberg time Th — (^f^i-i ■ Here 
fi(2?) is the phase volume density including spin degrees 
of freedom at the energy E, and / is the spacial dimen- 
sion. Throughout this paper, we study the two dimen- 
sional case / = 2. Two leads are assumed to have Ni and 
N2 channels, i.e. a± = 1, 2, • • • , JVi and a 2 = 1, 2, • • ■ , N 2 . 



The classical action of the orbit a is S a — J a p-dq, where 
q and p are the position and momentum variables. 

The stability amplitude is decoupled into two factors 
A a and A a . The first factor A a accounts for the sta- 
bility in the position and momentum space, and the sec- 
ond factor A Q originates from the spin dynamics. Both 
A a and A Q are uniquely determined, when the classical 
trajectory a governed by the microscopic Hamiltonian 
is specified. However their statistical behavior is inde- 
pendent of the details of the trajectory. As discussed in 
Refs l2lll22l and[H, the stability amplitude A a is related 
to the survival probability in the chaotic cavity. Postu- 
lating an ergodic motion in the position and momentum 
space, we obtain the following sum rule 

£ |-4 Q | 2 = f 'dTe-^/W (5) 

a.:a\— 

where T H /2N (N = Ni + N 2 ) can be regarded as the 
dwell time inside the cavity, so that the inverse is the 
escape rate. The escape rate is related to the position 
and momentum variables and is unrelated to the spin 
variables. Hence the dwell time should be Th/2N rather 
than Tjj/N. The spin matrix A Q is defined as 

A(i) = e #(i)^/ 2 e '«M^/2g#(t)^/2 (g) 

along a trajectory a, where <x = (o Xl a y , a z ) consists of 
the Pauli matrices. The time evolution of the Euler an- 
gles (ip(t), 0(t), 4>(t)) is microscopically determined by the 
Schrodinger equation 

ih^A(t) = HA(t), (7) 

where % is the effective Hamiltonian which describes the 
spin-orbit interaction. We assume that the spin dynam- 
ics is subdominant in the semiclassical limit. That is, 
the dynamics of the position and momentum variables 
are determined by the spacial Hamiltonian without spin 
degrees of freedom, while the spin is influenced by the 
momentum motion via the spin-orbit interaction. The ef- 
fective Hamiltonian in (J7J describes such a subdominant 
dynamics of the spin variables. A similar hierarchical 
structure has successfully been employed to analyze the 
GOE-GUE crossover regim o 22 i 28 ~ — : the resulting phys- 
ical quantities are in agreement with the corresponding 
RMT expressions. An RMT prediction for the paramet- 
ric spectral form factor in the GSE symmetry class was 
also reproduced in a similar way 3 ^. 

Since bumpy boundaries of a cavity induce a chaotic 
behavior in the position and momentum variables, the 
momentum effectively plays a role of a stochastic mag- 
netic field 3 -!. Then the effective Hamiltonian which de- 
scribes the time evolution of the spin can be written as 

U = lso h-(^o^, (8) 

where 7 so is the coupling constant and h = 
(h x (t),hy(t),h z (i)) is the effective stochastic magnetic 
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field. We assume that the classical spin undergoes a 
Brownian motion on the Bloch sphere^ due to the 
stochastic magnetic field satisfying 

((h a (t)h a ,(t'))) = 2V8 a>a ,8{t-t'), a,a' = x,y,z, (9) 

where the brackets ((...)) denote an average over the 
stochastic process of the magnetic field and T> is the 
diffusion constant. Then the probability density func- 
tion of the Euler angle P(ip,9,<fi) (with the measure 
sin 6di/jd8d(j)) obeys the Fokker-Planck equation 
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where C is the Laplace-Beltrami operator 
C = 
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This stochastic dynamics can exactly be analyzed by us- 
ing Wigner's D function^. Let us suppose that the ini- 
tial Euler angles at t = are u/ = (ip' , 9' , <f>'). Then the 
conditional probability to find the angles at oj = (-0, 9, (f>) 
after time t is 

°° j 3 •>; 4- 1 

j—0 m——j n——j 

x Dl^, 9, 4>){Dl in (i,', 9', 0')}* e~^ +1 ^ m . 

(12) 

Here an asterisk signifies a complex conjugate and 
Wigner's D function is defined as 
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in terms of the Jacobi polynomials Pjf' b \x). The index 
j is an integer or a half odd integer (j = 0, 1/2, 1, 3/2, • • ■ 
and m,n — — j, — j + 1, ■ ■ • The conditional probabil- 
ity g(aj;t\uj') satisfies a normalization condition 



dco g(ui; t\oo') — I d9 

lo Jo 



defy / dip sin 6* g(uj; t\uj') = 1. 



III. TRANSPORT IN THE GOE-GSE 
CROSSOVER REGIME 

A. Average conductance 

The average conductance G is written in terms of the 
transmission amplitude t ai a2 as 
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Here the transmission matrix t is a 2N% x 2N2 matrix 
which consists of the 2x2 blocks t ai ,a 2 - Then a semi- 
classical expression 
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follows from Here the brackets (...) mean an energy 
average, which eliminates the fluctuations of the physical 
quantities. If the difference between the actions S a and 
Sj is sufficiently large, the exponential term e^^"^ 3 ^ 
rapidly oscillates in the semiclassical limit h — > 0, which 
eventually vanishes after averaging. Hence, in order to 
give a finite contribution, the trajectories a and 7 are mu- 
tually almost the same. Then the identical trajectories 
a = 7 yield the first order approximation, which is re- 
ferred to as "the diagonal approximation"^. These mu- 
tually identical trajectory pairs yield the following con- 
tribution 

(TrCttt))! = A^^Tr(A Q At)|^| 2 
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Here we used the sum rule ([5]) for the stability amplitude 
A a . In this calculation, a product of the spin matrices is 
reduced to the identity matrix, and the trace yields a fac- 
tor 2. The diagonal approximation does not discriminate 
the symmetry classes. 

The second order approximation comes from the 
Richter-Sieber (RS) pair o 26 i 27 drawn in Fig.l. In the RS 
pair, two trajectories come close to each other in the en- 
counter region, and go in the opposite directions on one 
loop. We can symbolically write RS pairs (See Fig.l) as 

a P\EL 2 EL 3 , 
7 : LiEL 2 EL 3 , 
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where E implies one of the two trajectory segments in 
the encounter region where two loops are connected and 
E implies the time reverse of E. The loops are denoted 
as Li, L2 and L3, respectively, and Lj {j — 1, 2, 3) is the 
time reverse of Lj. Using these notations, we can write 
the RS pair contribution as 
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where AS* is the action difference S a — 5 7 . The spin 
matrices A Q and A 7 are factored into the loop and en- 
counter parts as 

A Q 
A, 



A Ll A £ A L2 (A B )- 1 A L 
A Ll A E (A L2 )- x (A E )-'A 



L 3 - 



(19) 
(20) 



Along the trajectories the non-Abelian nature of the spin 
operators must be taken into account. A time reversal 
operation of a spin matrix is realized by a matrix inver- 
sion. 




FIG. 1: The Richter-Sieber (RS) pair. The solid and dashed 
curves are respectively a and 7 orbits in the text. 

We divide the whole time elapsed on a trajectory into 
the loop and encounter parts, i.e., Ti,T 2 ,T 3 for L\, L 2 , 
L3, respectively, and t enc for E. It should be noted that 
the existence of the encounter affects the survival prob- 
ability e -( 2N / T H)T m Eq> @21,22. if Qne of the tw0 Qr _ 

bit segments in the encounter is inside the cavity, the 
other segment must also remain inside. Hence the sur- 



vival probability is modified into e 



-(2N/T H )(T-t snc ) 



In the encounter region, the classical actions of the two 

trajectories are slightly different. The action difference 

can be estimated by using the coordinates (s, u) along the 

stable and unstable manifolds within the ranges s,u € 

[— c, c] 21 ' 22 . The time duration t enc inside the encounter 

region is related to the Lyapunov exponent A as t cnc ~ 
1 2 

v In t^t , and the action difference is expressed as 
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(21) 



The number density of encounters in a trajectory with 
an elapsed time T = Ti + T 2 + T 3 + 2i cnc is evaluated 
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Ti,T 2 >0 
' T 1 +T 2 <T-2t on 

by taking account of T u T 2 > and T 3 = T - T x - T 2 

2^nnn ^ 0- 



On the other hand, the spin diffusion term is calculated 



as 
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Using the above formulas, we find that the RS contribu- 
tion is 
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The last line of the above formula is obtained by the fol- 
lowing criterion: after expanding the formula in i e nc and 
integrating each term over (s, u), any term dependent on 
t Gnc vanishes in the semiclassical limit, and a finite con- 
tribution comes only from the terms independent of t enc . 

The third order contribution to the average conduc- 
tance comes from the diagrams drawn in Fig. 2. In Ap- 
pendix A, the spin diffusion terms are listed. Using these 
results, we arrive at an expression of the average conduc- 
tance 
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(25) 



where £ is the crossover parameter defined as £ = 
Js Q VTh/N. We can easily check that (f25|) reproduces 
the 1/N expansions of the GOE and GSE formulas ([T]) 
by taking the limits £ — > and £ — > oo, respectively. 
It follows from this result that only one parameter £ is 
necessary to describe the GOE-GSE crossover. 
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FIG. 2: The diagrams of the third order contributing to the conductance 



B. Conductance variance 

The conductance variance VarG is written in terms of 
the transmission matrix t as 



VarG 



= ({Tr(ttt)} 2 ) - (Tr(ttt)) 2 . (26) 



Let us first adopt the diagonal approximation, for which 
we have two types of contributions. One is given by set- 
ting a — f3 and 7 = S, where ai,a 2 ,ci and c 2 are all 
independent. The other choice is a = 5, 7 = /3, and 
di = Ci , C12 = C2 . The numbers of possible channel com- 
binations are NfN$ and NiN 2 , respectively. These con- 
tributions are summed up to yield 



The semiclassical expression of the first term is 
,2\ 4 
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Here the spin diffusion term is obtain the diagonal contribution 

«|Tr(A a At)| 2 )) , 3 v _ 4N?NI N lN 



{Tr(ttt)}^) = + i^L 1 + . (30) 



= 1 + 3e~ 2j * a ' D ( To ' +T ~<\ (29) To go beyond the diagram approximation, we note that 

the next order diagrams are classified into rf-families and 
where T a and T 7 are the times elapsed on the trajectories x-families as shown in Fig. 3^2: (i-quadruplets are drawn 
a and 7, respectively. Using this expression, one can in the diagrams (i)-(vii), while ^-quadruplets in (viii). 




FIG. 3: The diagrams contributing to the conductance variance 



Let us write the next order term as 
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Here the coefficients di, d 2 and x\ are obtained from the 
families of quadruplets. The coefficient d\ comes from the 
diagram (i): quadruplets consisting of one diagonal pair 
and one RS pair. On the other hand, many diagrams 
have to be taken into account to calculate d 2l i.e., 1): 
quadruplets consisting of one diagonal pair and one pair 
contributing to the 0(1/N) term in the expansion ([25]) 
of the average conductance, 2): two RS pairs 3): the 
diagrams (ii)-(vii) shown in Fig. 3. The coefficient x\ is 
calculated from the diagram (viii) in Fig. 3. Considering 



the contributions from these diagrams, we obtain 
12 
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Then we find the expression of the conductance variance 
for iVi > 1 and N 2 > I as 



VarG _ NjNl f 6 \ (1 



One can check that the GOE and GSE limits agree with 
©■ 
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C. Shot noise 

The shot noise P is written in the form 

— = (Tr(tt t -tt t tt t )). 
"o 

The second term is semiclassically expressed as 
4 



(32) 
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x (^A Q A}A 1 A* s er^ s "- s f+ s -'- s ^ 

x ((Tr(A a AtA 7 Al))}). (33) 

The diagonal contribution consists of two terms. One 
term has a = ft and 7 = 5 where we need to set a\ = c%. 
The other term has a = 6, /3 = 7, and <Z2 = C2. We sum 
up these two terms and obtain 
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(34) 



It follows that the diagonal contributions to (Tr(tt'l')) and 
(Tr(tttttt)) are both 2^^-, and mutually cancel. 

Hence the RS pair contribution to Tr(ttt), and the 
contribution to Tr(tt^tt'l') from the quadruplets drawn 
in Fig. 3 (i) and (viii) have to be calculated. Moreover 
we take account of the additional diagrams shown in 
Fig. l 21 i 22 . Summing up these contributions, we finally 
obtain the shot noise in the crossover regime as 
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(35) 



Let us denote the O(l) terms of G and P by 5G and SP, 
respectively. It can be seen from Eqs. (|2"5)l and (|3"5|) that 
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which is a universal relation established in Refll_j 
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FIG. 4: The diagrams contributing to the shot noise 
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IV. COMPARISON WITH A RANDOM 
MATRIX MODEL 

In this section, a random matrix model on the GOE- 
GSE crossover is numerically analyzed and the results are 



compared with the semiclassical formulas. In the theory 
of random matrices, a time-reversal invariant quantum 
Hamiltonian with spin 1/2 is simulated by a self-dual 



8 



real quaternion random matrix^!. A real quaternion q is 
a linear combination 



q = q e + q\&\ + qie-i + 9363, 



(37) 



where qj are real numbers and called the j-th component 
of q. The bases eo, e\, can be represented by 2 x 2 

matrices as 
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(38) 



so that qoeo is equated with a real number go- The dual 
quaternion of q is defined as 



q = qoeo ~ 9i ei - <?2e 2 - q^z- 



(39) 



When an m x n real quaternion matrix Q has the (j, /) 
element qji, we define that the n x m dual matrix Q 
has the element gjy. If a square real quaternion 

matrix satisfies Q — Q, then Q is called a self-dual real 
quaternion matrix. 

The parametric motion of a self-dual real quaternion 
random matrix is realized in the framework of Dyson's 
matrix Brownian motion models. It is postulated that 
the probability density function (p.d.f.) of an M x M 
self-dual real quaternion matrix H is 

P(H; t\R) dH oc exp |-2 ^f ~ | dH (40) 



1-e- 



with 
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M 3 



dH = HdH jj Hl[dH^. 



(41) 



j'=i 



Here ffj* is the fe-th component of ffjj. We are inter- 
ested in the parametric motion of the matrix H depend- 
ing on the fictitious time parameter r. 

At the initial time r = 0, this p.d.f. is reduced to 



P(H; 0\R) dH = 8{H - R) dH, 



(42) 



so that the self-dual real quaternion matrix R gives the 
initial condition of the parametric motion. On the other 
hand, in the limit r — > 00, 

P(H; oo\R) dH cx exp (-2Tr# 2 ) dH, (43) 

which is the p.d.f. of the GSE. 

Let us suppose that the elements of the initial matrix 
R have only the 0-th components (R is then a real sym- 
metric matrix) and that the p.d.f. of R is 

Pgoe (R)dR cx exp ( - ^Tri? 2 J dR (44) 



with 
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(45) 



which is the p.d.f. of the GOE. Then we can calculate 
the p.d.f. of H at a fictitious time r as 

P{H)dH = jy dR P{H;t\R)P goe {R)^ dH 
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which describes the crossover from the GOE (at r = 0) 
to GSE (at r = 00). The components of the elements 
Hji (j < I) are independently distributed according to 
Gaussian density functions. 

If an M x M real quaternion matrix Z satisfies 



ZZ = ZZ 



(47) 



(Im is the M x M identity matrix), then Z is called a 
symplectic matrix. If the elements of a symplectic matrix 
U only have the 0-th elements, then U is a real orthogo- 
nal matrix. It is known that the measure dH is invariant 
under the symplectic transformation H <— ¥ ZHZ and the 
measure dS is invariant under the orthogonal transforma- 
tion S H> U T SU (U T is the transpose of U). It follows 
that the p.d.f. P(H)dH in (|4"6"| is invariant under the 
orthogonal transformation H H> U T HU. 

We go back to the problem of a chaotic cavity with M 
bound states to which two leads with N\ and N2 propa- 
gating modes are attached(M > N — Ni + N 2 ). We are 
interested in the limit M — > 00. Let us suppose that the 
M x M matrix H describing the scattering in the cavity 
is a random matrix distributed according to the crossover 
p.d.f. in (|4*o| . Then the N x N scattering matrix S is2 



S = I 



N 



i2irW T {H - E F - iitWW T )- 1 W, (48) 



where Ep is the Fermi energy, and the elements of an 
M x N real matrix W are the coupling constants between 
the cavity and the leads. 

Assuming that the tunnel probability of the leads are 1, 
we can see that the eigenvalues of W T W are all M/(p7r 2 ), 
where p is the eigenvalue density of H at the Fermi en- 
ergy. Then a singular value decomposition 



W = UDV 



(49) 
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holds, where U and V are M x M and N x N real orthog- 
onal matrices, respectively, and D is an M x N matrix 



with 



D 



W 




In 




(50) 



(51) 



Here O is an (M — N) x N matrix consisting of zero ele- 
ments. When the Fermi energy Ep is set to zero (so that 
p = y/2M/n), the scattering matrix S can be rewritten 
as 



S = I N + W2MW 1 yH 
where 



1 M ~ - T 
-WW T 

2 



H = (UV) T H(UV) 



W, (52) 



(53) 



with 



V 



V o 

O Im-n 



(54) 



Since UV is a real orthogonal matrix, the M x M matrix 
H is also distributed according to the p.d.f. in (14*51) (with 
H replaced by H). 

Thus the scattering matrix S can numerically be gen- 
erated by using the Gaussian p.d.f. (|4"o"[) of H and the 
relation (1521) . Replacing the quaternion elements of S by 
the 2x2 matrix representations (|38|) . we obtain a 2Nx2N 
matrix S. It is written in terms of the 2Ni x 2N2 trans- 
mission matrix t as 



S 



r tt 
t r' 



(55) 



where r and r' are the reflection matrices. Then the 
average conductance G can be evaluated as 



G_ 

Go 



/2(N 1 +N 2 ) 2N 2 

(Tr(ttt)) = ( Yl El^l 

\j=2N 2 + l 1=1 



(56) 



where the brackets (...) denote an average over the p.d.f. 
(fifl)) . The conductance variance VarG and the shot noise 
P can similarly be written as 

^ = ({Tr(ttt)} 2 \-(Tr(ttt))2 (57) 



and 



P 
Pn 



(Tr(ttt) - Tr^tttt 1 )) 



(58) 



The remaining task is to find a relation between 
the semiclassical parameter £ and the fictitious time t. 



Pandey analyzed hierarchical relations among the eigen- 
value correlation functions of random matrices and eval- 
uated the form factor K{k\r) (the Fourier transform of 
the scaled two eigenvalue correlation function). For the 
random matrices obeying the crossover p.d.f. in (|46[) . he 
derived a relation^ 

K(k; t) = K(k; 00) + {K(k; 0) - K(k; 00)} e" 4 * 1 '' 2 "' 

(59) 

with k I 0. Here K(k; 0) and K(k; 00) are the form fac- 
tors of the GOE and GSE random matrices, respectively. 
They are known to be^i 



K(k;0) = 2k, K(k; 00) 



so that 



^(fc;r) = -(l + 3e 



-8Mrk\ 



k I 0, 



(60) 



(61) 



On the other hand, Nagao and Saito semiclassically an- 
alyzed the form factor of a chaotic system with a weak 
spin-orbit interaction. They obtained a small k > ex- 
pansion up to the second order— 

K(k;a) = ^(1+3 e - akT »)+^ {l + (3akT H - 9) e ~ akTH } 

with a — 7s X>. Comparing (f6T1) and ([(52^) . we arrive at a 
relation 



8Mt = aT H . 



(63) 



Therefore the semiclassical parameter £ = aTjj/N is as- 
sociated with the random matrix parameter r as 



8M r 
~N~'' 



(64) 



8.15 
8.1 
8.05 



7.95 



G/G„ 



7.9 
7.85 

7.8 
7.75 

7.7 
7.65 



0.02 



0.04 



0.06 



0.08 



0.1 



FIG. 5: A comparison of the semiclassical result (curve) and 
numerically calculated random matrix results (errorbars) for 
G/G Q . 
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In Fig. 5, numerical calculations of G/Go at various val- 
ues of t are compared with the corresponding semiclassi- 
cal predictions (f25j) with £ = 8Mt/N. In the numerical 
calculations, we set M = 200, N x = 20 and N 2 = 5. The 
crrorbars arc introduced in order to estimate the statisti- 
cal errors due to the fact that the averages are calculated 
over only 300 samples. Note that the semiclassical for- 
mulas are truncated and hence are valid only for large 
Ni and N 2 . Nevertheless we can see a fairly reasonable 
agreement with the numerical results. 

Similar plots for VarG/G 2 , and P/Po are also shown, 
respectively, in Figs. 6 and 7. The semiclassical curves 
are drawn by using eqs. (|3"Tj) and (|35[) . Since we again 
find reasonable agreements, it can be conjectured that 
there is an equivalence between the semiclassical method 
and random matrix theory. 

In the case of the GOE-GUE crossover, the correspond- 
ing random matrix model can analytically be treated^ 
and the results can be compared with the semiclassical 
formulas 22 . It seems possible to apply similar techniques 
to the GOE-GSE crossover. For example, the diagram- 
matic perturbation theory is able to give the leading 
terms of the transport properties^—. It would be inter- 
esting to compare the semiclassical formulas with such 
analytical results and confirm the equivalence mentioned 
above. 



0.25 



0.15 - 



VarG/G, 




0.05 



0.04 X 006 



FIG. 6: A comparison of the semiclassical result (curve) and 
numerically calculated random matrix results (errorbars) for 
VarG/Gg. 



V. SUMMARY 

We studied a chaotic quantum transport of an electron 
with spin-orbit interaction in a cavity. Our approach is 
based on the semiclassical theory. The key ingredient 
of the theory is the universal statistics of the stability 
amplitudes. The electron diffusion in the position and 
momentum space is related to the escape rate, and the 



1.4 
1.38 
1.36 
1.34 
1.32 

1.3 
1.28 
1.26 
1.24 
1.22 




0.02 



0.04 



0.06 



0.08 



0.1 



FIG. 7: A comparison of the semiclassical result (curve) and 
numerically calculated random matrix results (errorbars) for 
P/Po- 



spin diffuses on the Bloch sphere due to the spin-orbit in- 
teraction, where the momentum variable plays the role of 
a stochastic magnetic field. Consequently the crossover 
parameter depends on the diffusion constant of the spin. 
The spin diffusion terms appear in a non-Abelian way 
along the classical trajectories. The spins along the tra- 
jectories interfere with each other, resulting in the change 
of the total spin. For instance, Eq. (|2"31 has both singlet 
and triplet contributions. These kinds of interference ef- 
fects seem to play a crucial role. 

In our calculation of the physical quantities such as 
the average conductance, only the first several terms of 
the resulting expansions were worked out. Such trun- 
cated results are valid only when the channel numbers 
are large. In order to obtain the full expansions, a more 
systematic calculation of the spin diffusion terms would 
be necessary. Although our expansions are truncated, 
they are still useful in the experimental point of view, 
because the channel numbers can be very large. More- 
over the GOE-GSE crossover can be realized when the 
spin-orbit interaction is controlled by applying an elec- 
tric field. Therefore we believe that an experimental test 
of our theory is in principle possible. 

In addition to the transport properties analyzed in this 
paper, the shot-noise variance is also an important quan- 
tity, which can be treated in an RMT framework^ -. It 
seems possible to apply the semiclassical method to the 
shot-noise variance. More ambitiously, one might be able 
to calculate arbitrary order cumulants of the conductance 
and shot noise in a semiclassical framework, as the RMT 
approach has already made a progress in that directions^. 

Our theory is valid in the case that the dwell time 
of the electron is much larger than the Ehrenfest time. 
When the Ehrenfest time is relatively large, the resulting 
corrections should be considered^ 2 -. In addition, the spin 
diffusion mechanism might depend on the specific form 
of the spin-orbit couplin g 15 ' 20 . These problems are also 
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interesting in experimentally realizable situations. 
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Appendix A: Spin diffusion terms 



In this Appendix, we present the spin diffusion terms contributing to the average conductance, conductance variance 
and shot noise. Here Tj and tj are the times elapsed on the loop Lj and encounter Ej, respectively, and a = 7g £>. 
These spin diffusion terms were evaluated by using Mathematica. 
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1. Average conductance 



The spin diffusion terms contributing to the average conductance (corresponding to the diagrams shown in Fig. 2) 
are fisted below. 

(i) : ^Tr | A Ll A El A L2 A E2 A Ls A El A Li A E2 A L& (A Ll A Bl A i4 A B2 A i3 A Bl A La A B2 A L5 ) f }^ 
= ((Tr{A L2 A £2 A i3 A £l A i4 (A L2 )- 1 (A £l )- 1 (A L3 )- 1 (A £2 )- 1 (A i4 )- 1 })) 

_ 1 _|_ 3 e -a(2ti+2t 2 +2T 2 +2T 3 ) _|_ 3 g-a(2T 2 +2T 4 ) _|_ 3 e -a(2t 1 +2t 2 +2T 3 +2T 4 ) _ £ e -a{2t 1 +2t 2 +2T 2 +2T 3 +2T 4 ) _ 



(ii) : ({Tr [A Li A Ei A L2 A E2 A L3 (A Ei )- 1 ^l 4 A E2 A L5 {A Ll A £l (A L3 )- 1 (A iS2 )- 1 (A i4 )- 1 A Bl A L2 A B2 A L5 } 1 
= ((Tr{A i2 A B2 A L3 (A £l )- 1 A L4 (A i2 )- 1 (A £l )- 1 A L4 A B2 A L3 })) 



(hi) 



I - l e -a(2t 2 +2T 2 +2T 3 ) _ | e -a(2t 1 +2T 2 +2T 4 ) _|_ | e -a(2t 1 +2t 2 +2T 3 +2T 4 ) _|_ 3 e -a(2t 1 +2t 2 +2T 2 +2T 3 +2T 4 ) _ 

Tr [A Ll A El A L2 A E2 A L3 A El A La ( A B2 ) ~ 1 A Ls 

x {A il A £l A i4 (A £2 )- 1 (A L2 )- 1 (A Bl )- 1 (A i3 )- 1 (A £2 )- 1 A i5 } t ] )) 
((Tr{A L2 A E2 A L3 A iSl A L4 A L3 AB 1 A L2 A £2 (A L4 )- 1 })) 

1 _ | e -a(2t 1 +271,+2T 4 ) _|_ | e -a(2t 1 +2t 2 +2T 2 +2T 3 ) _ 3 g -a(2t 2 +2T 2 +2T 4 ) _|_ 3g-a(2ti+2t 2 +2T 2 +2T 3 +2T 4 ) _ 



(iv) 



(v) 



Tr [ A L 1 A El A L2 ( A El ) - 1 A La A E2 A i4 ( A E . 2 ) ~ 1 A Ls 

x {A Ll A £l (A L2 )- 1 (A £l )- 1 A i3 A B2 (A L4 )- 1 (A £2 )- 1 A L5 } t 
-((Tr{A i2 (A £l )- 1 A i3 A B2 A i4 A L4 (A £2 )- 1 (A i3 )- 1 A Bl A i2 })) 

__ I _ 3 p -a2T 2 , 9 -a(2T 2 +2T 4 ) _ 3 p -a2T 4 
2 2 2 2 

Tr A Ll A El A L2 A E2 A L;i (Ae^Al, (A El )" 1 A L , 

x {A il A Bl (A i4 )- 1 A B2 A L3 (A B2 )- 1 (A L2 )- 1 (A £l )- 1 A i5 } f 
= ((Tr{A L2 A B2 A L3 (A £2 )- 1 A L4 A i2 A B2 (A i3 )- 1 (A B2 )- 1 A L4 })) 

_ 1 _ 3 -a2T 3 i 3 -o(2T 2 +2T 4 ) , 3 -a(2T 2 +2T 3 +2T 4 ) 
22 ^2 ^2 



(vi) : ^Tr^.A^A^A^A^A^A^IA^.A^A^A^A^A^AL,} 1, 
= ((Tr{A i2 A £l A L3 (A L2 )- 1 (A £l )- 1 (A L3 )- 1 })) 

— I _|_ 3 e -a(2t!+2T 2 ) _|_ 3 e -a(2t!+2T 3 ) _|_ 3 g-a(2T 2 +2T 3 ) _ 3 e -a(2i 1 +2T 2 +2T 3 ) _ 



(vii) : ((Tr [A il A £l A i2 (A Bl )- 1 A i3 A Bl A L4 {A Ll A iJl (A L2 )- 1 (A iSl )- 1 (A i3 )- 1 A Bl A L4 } 1 

= (( Tr { A L2 ( A El ) - 1 A L3 A L;i A El A L2 } )) 

_ 1 _ 3 p -a2T 2 i 9 -o(2T 2 +2T 3 ) _ 3 -a2T 3 
2 2 ' 2 2 

(viii) : ^Tr [a Li A Bi A L2 A Ei A L3 (A Ei )- 1 A L4 {A Ll A El A L3 (A Bl )- 1 (A L2 )- 1 (A £l )- 1 A L4 } t 

= ((Tr{A L2 A £l A L3 A L2 A £l (A i3 )- 1 })) 

_ I i 3 -o(2ti+2T 2 ) , 3 -a(2ti+2T 2 +2T 3 ) _ 3 -a2T 3 



(ix) 



Tr 



A Ll A Bl A L2 (A iJl )- 1 A L3 (A iJl )- 1 A L4 {A Ll A Bl (A L3 )- 1 A iSl A L2 (A iSl )- 1 A L4 } t 



<(Tr{A i2 (A £l )- 1 A i3 (A L2 )- 1 (A £l )- 1 A i3 })) 

1 _ 3„-a2T 2 , 3 -a(2ti+2T 3 ) , 3 -a(2ti+2T 2 +2T 3 ) 

2 2 c "T 2 "I" 2 



13 



2. Conductance variance 

The spin diffusion terms contributing to the conductance variance (corresponding to the diagrams shown in Fig. 3) 
are listed below. 



(i) : ((Tr A Ll A El A L2 (AeJ-'A^ {A^A^A^r 1 ^)-^}^ Tr{A L4 (A i4 )t} 

= ((2Tr{(A i2 ) 2 })) = -2 + 6e- a2T2 . 

(ii) : ((Tr {a Li A Ei A L2 A E2 A L3 (A Ll A El A l , A B2 A L3 ) t } Tr {a L4 A El A l , A E2 A L(i (A Li A El A L2 A E2 A L6 ) f } 
|Tr { A L2 ( A i5 ) t } | 2 )) = 1 + 3e- a ^+ 2T ^ . 



(hi) : ((Tr j Al 1 A El Al 2 A E2 Al 3 A El Al 4 (A Ll A El A L J f } Tr {a Ls A E2 A L(i (A^ 5 A E2 A Es A El A E2 A E2 Al 6 )^ | ^ 



= (( |Tr(A L2 A £2 A L3 A £l )| 2 )) = 1 + 3e -( 2 ™+ 2 ^+ 2 ^). 



(iv) : ((Tr {A Ll A El A L2 A El A L3 (A Ll A El A L3 ) t }Tr{A i4 A £l A L5 (A i4 A £l A i2 A Bl A i5 )t} 
Vr(A L2 A £l )| 2 \\ = l + 3e- a ( 2tl + 2T2 ). 



(v) : ((Tr [A Ll A £l A i2 A B2 A i3 {A Ll A^A^^A^A^} 1 

A i4 (A B2 )- 1 A i5 (A £l )- 1 A L6 {A L4 (A £2 )- 1 (A i2 )- 1 (A Bl )- 1 A L6 } t 
{Tr(A i2 A L5 )} 2 \\ = 1 + 3 e - a < 2T2+2T =). 



x Tr 



(vi) : ((TrjA^A^A^A^A^A^A^ (A Ll A^ AlJ 1 } 

x Tr [a^Ae^Al, {A i5 (A £2 )- 1 (A i2 )- 1 (A Bl )- 1 (A L3 )- 1 (A i?2 )- 1 A i6 } t ' 

'{T^A^A^A^)}^ = 1 +3e -a(2t 1+ 2t 2+ 2T 2+ 2T 3 )_ 



(vii): ((Tr{A Ll A £l A i2 A £l A i3 (A Ll A Bl A L3 ) t }Tr[A L4 (A iSl )- l A L5 {A i4 (A Bl )- l (A i2 )- l (A Bl )- l A L5 } 1 
= (({Tr(A i2 A Bl )} 2 \\ = l + 3e- a ( 2tl+2T2 ). 



(viii) : ((Tr{A Ll A Bl A L2 (A Ll A El A L J f } Tr {a l , 3 A El A Li (A Ls A El AlJ*})) 



Tr(A L2 A^ 4 ) 



= 1 + 3 e - a ( 2T 2+2T 4 )^ 



3. Shot noise 

The spin diffusion terms contributing to the shot noise (corresponding to the diagrams shown in Fig. 4) are listed 
below. 



(i) : ((Tr |a Li A £i Az, 2 (A L; , A^ A L2 ) f A^A^A^A^A^A^) 1 A U 

x {A Ll A Bl A L4 A £2 (A L5 )- 1 (A £2 )- 1 A L6 } t 
= ((Tr{(A L5 ) 2 })) = -l + 3 e - 2 ^. 



14 



Tr 



A Ll A El A L2 {A i6 A Ea (Ai 3 )- 1 AB 1 A L2 } t A Ls A £2 A L4 (A Bl )- 1 A L3 (A B2 )- 1 A L6 
{A^A^A^)- 1 ^)-^} 1 ])) 
Tr{(A^A i3 A i4 ) 2 })) = -l + 3 e - 2a ( tl+ ^). 



Tr 



A^A^A^ (A L3 A El A L2 ) f Al.Ae.Al^AeJ^Al, {A Ll A i3l (A L4 )- 1 (A jEl )- 1 A L5 } t 



= ((Tr{(A L4 ) 2 })) = -l + 3 e - 2 ^. 



